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THE MOLECULAR GAS CONTENT OF Z < 0.1 RADIO GALAXIES: 
LINKING THE AGN ACCRETION MODE TO HOST GALAXY PROPERTIES 

V. Smolcic^'^'^ and D. a. Riechers^"'' 
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ABSTRACT 

One of the main achievements in modern cosmology is the so-called 'unified model', which success- 
fully describes most classes of active galactic nuclei (AGN) within a single physical scheme. However, 
there is a particular class of radio-luminous AGN that presently cannot be explained within this 
framework - the 'low-excitation' radio AGN (LERAGN). Recently, a scenario has been put forward 
which predicts that LERAGN, and their regular 'high-excitation' radio AGN (HERAGN) counter- 
parts represent different (red sequence vs. green valley) phases of galaxy evolution. These different 
evolutionary states are also expected to be refiected in their host galaxy properties, in particular their 
cold gas content. To test this, here we present C0(1— >0) observations toward a sample of 11 of these 
systems conducted with CARMA. Combining our observations with literature data, we derive molec- 
ular gas masses (or upper limits) for a complete, representative, sample of 21 z < 0.1 radio AGN. 
Our results yield that HERAGN on average have a factor of ~ 7 higher gas masses than LERAGN. 
We also infer younger stellar ages, lower stellar, halo, and central supermassive black masses, as well 
as higher black hole accretion efficiencies in HERAGN relative to LERAGN. These findings support 
the idea that high- and low-excitation radio AGN form two physically distinct populations of galaxies 
that reflect different stages of massive galaxy build-up. 

Subject headings: galaxies: fundamental parameters - galaxies: active, evolution - cosmology: obser- 
vations - radio continuum: galaxies 



1. INTRODUCTION 

Over the past two decades a standard model of AGN 
has emerged. In this 'unified' model efficient disk accre- 
tion of cold matter on the central supermassive black 
hole (BH) provides the radiation field that photoion- 
izes emission-line regions. However, there is a certain 
fraction of AGN identified by radio observations that 
poses a challenge to the unified model, the so-called low- 
excitation radio AGN (hereafter: LERAGN). The main 
difference between high- excitation radio AGN (HER- 
AGN) and these LERAGN is that the latter do not ex- 
hibit strong emission lines in their optical spectra (Jack- 
son & Rawlings 1997; Evans et al. 2006). 

Recently, Hardcastle et al. (2006) have suggested that 
high- and low-excitation radio AGN may represent a 
principal separator between populations fundamentally 
different in their black hole accretion mechanisms (see 
also Evans et al. 2006; Allen et al. 2006; Kewley et al. 
2006). They developed a model in which central super- 
massive black holes of HERAGN accrete in a standard 
(radiatively efficient) way from the cold phase of the in- 
tragalactic medium (IGM), while those of LERAGN are 
powered in a radiatively inefficient manner by Bondi ac- 
cretion of the hot IGM. Smolcic (2009) showed that low- 
and high excitation radio AGN exhibit not only systemic 
differences in their black hole masses and accretion rate 
properties, but also in their host galaxy properties, such 
as stellar masses and stellar populations. This is consis- 
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tent with these two classes of radio AGN dividing in a 
stellar mass vs. color plane in such a way that LER- 
AGN occupy the red sequence and HERAGN inhabit 
the so called "green valley" , a sparsely populated region 
between the blue-cloud and the red-sequence (Smolcic 
2009). 

The stellar mass vs. color plane can be interpreted 
as a time-sequence for galaxy evolution. Galaxies 
are thought to evolve from an initial star-formation- 
dominated state with blue optical colors into the most 
massive red-and-dead galaxies through a transition phase 
reflected in the green valley (Bell et al. 2004a, 2004b; 
Borch et al. 2006; Faber et al. 2007; Brown et 
al. 2007). In recent years it has been suggested 
that radio outflows from AGN lik ely play a crucial 
role in this massive galaxy build-up (Croton et al."200 ' 
Bower et al. 200 6; SUacki fc Springel 2006; Siiacki et al 
2007t iFanidaki s et al.ll2010t ). In this context the radio- 
AGN feedback (often called the "radio" or "mainte- 
nance" mode), which is thought to limit stellar mass 
growth in already massive galaxies, is expected to occur 
only in LERAGN (Smolcic 2009). 

Furthermore, it has been shown that the cosmic evolu- 
tion of the space density of various types of radio AGN 
is significantly different (e.g. Peacock et al. 1985, Willott 
et al. 2001; Smolcic et al. 2009). Based on a study of the 
evolution of the radio AGN luminosity function out to 
z = 1.3, Smolcic et al. (2009) have shown that the comov- 
ing space density of low-luminosity radio AGN (predom- 
inantly LERAGN) only modestly declines since z — 1.3, 
while that of powerful AGN (predominantly HERAGN) 
dramatically diminishes over the same cosmic time inter- 
val. This suggests that LERAGN and HERAGN not only 
represent physically distinct galaxy populations, but also 
populations in different stages of massive galaxy build- 
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up. If this is the case, the molecular gas masses and 
fractions in low- and high- excitation radio AGN are ex- 
pected to directly reflect this trend. 

We here investigate this idea by observing 
C0(J=1— J>0) emission of a carefully selected, rep- 
resentative sample of nearby {z < 0.1) HE- and 
LERAGN with CARMA. We adopt a ACDM cosmology 
with H„ = 70, = 0.3, = 0.7. 

2. DATA 
2.1. Sample 

We here utilize a sample of 21 Type 2 AGN at z < 0.1 
that have been obse rved in X - rays (with Chandra or 
XMM-Newton) bv E vans et"all (|2006[ ). 18 out of the 21 
AGN have been drawn from the 3CRR survey, adding 
3 more sourc es (3C 403, 3C405 and Gen A) for com- 
pleteness (see lEvans et al.l 120061 for details) . The sam- 
ple properties are summarized in tab:physprops (see 
also Tab. 1 in Evans et al. 2006). We separate our 
AGN into LERAGN (i.e. LINERs) and HERAGN (i.e. 
Seyferts) using standard diagnostic tools based on opti- 
cal emission line flux ratios where possible (see Fig. [T] and 
tab:physpro ps ; [K auffmann et al.' '2003'; 'Kewle v et al.l 
[2001, 200(| ISmolei c 2009; Buttiglionc ct al. 2009). For 
this we make use of the emission line fluxes extracted 
from high resolution spectr oscopy of 3GR sources pre- 
sented in iButtiglione efaD (|2009l . 120101 I2011L see also 
Tab. 1 in iSmolcig |2OO0^ In cases where the relevant 
emission line fluxes are not available, we make use of 
the galaxy type information available in the NASA Ex- 
tragalactic Database0 to separate the sources into LE- 
and HE-RAGN. The sample contains 9 HERAGN and 
12 LERAGN. 

2.2. CO(l^O) observations and data reduction 

At the time of observations, 8 out of the 21 Type 2 
AGN in our sample have already been detected in 
GO(l-^O). Thus, we observed the GO(l-)-O) transition 
line toward the remaining 13 AGN using the the CARMA 
(Combined Array for Research in Millimeter-wave As- 
tronomy) Interferometer. Observations were performed 
during Summer/2009 and Spring/2010 for about 4 to 
15 hours per source (tabiobs ). All targets were ob- 
served under good to excellent weather conditions at 
3 mm with 15 antennas (corresponding to 105 baselines) 
in the two most compact, E and D configurations (2009 
and 2010, respectively). Data on two objects (3C 388 
and 3C 405) had to be discarded due to technical prob- 
lems, and are excluded in the following. The receivers 
were tuned to the redshifted C0(J=1— J-O) line frequen- 
cies (i/i-est=115. 2712 GHz; see tab:obs for exact observ- 
ing frequencies), centering them in the upper sideband. 
Three bands with 15 channels of 31.25 MHz width each 
were utilized. The bands were overlapped by 2 channels 
to improve calibration of the correlated dataset, leading 
to an effective bandwidth of 1281.25 MHz (3500 km s'^) 
per sideband. Phase calibration was performed by ob- 
serving bright nearby radio quasars every 15 minutes. 
Bandpass calibration was performed once per track on 
bright quasars. Fluxes were bootstrapped relative to 
planets, or monitored radio quasars if no planet was 
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available. The total calibration is estimated to be ac- 
curate to 15%. Data reduction was performed using the 
MIRIAD package. The C0(1— >0) spectra are shown in 
Fig. [2 

3. RESULTS 
3.1. CO Data 

CO(l^O) has been detected in 4 (3C33, 3C321, 3C403, 
and NGC6109) out of the 11 galaxies in our CARMA-CO 
sample (see Fig. [2]). To parameterize the emission lines 
detected in these four galaxies we fit Gaussian profiles to 
the line and underlying continuum emission (see tabxo 
for line/continuum properties). 

Two of our four CO detected sources (3C321 and 
3C403) have recently been detected in the CO(l^O) 
transition by Flaquer et al. (2010) using the IRAM 30 m 
telescope. The line parameters reported by Flaquer et 
al. are in good agreement with ours. 

3(7 upper limits for C0(1— !>0) non-detections are de- 
termined by assuming a line width of 300 km s~^, corre- 
sponding to the average width of the detected lines. We 
further complement tab:co with data from literature for 
the 8 sources with already existing C0(1— >0) detections, 
and the 2 sources (3C 388 and 3C405) that had to be 
excluded from our sample. 

3.2. Ancillary Data 

We summarize the physical properties of the 21 sources 
in our low- and high-excitation radio AGN sample in 
tabiphysprops . We adopt the 178 MHz luminosities, ac- 
cretion efficiencies, and black hole masses from Evans et 
al. (2006). The stellar ages of our sources, taken from 
Buttiglione et al. (2009), were derived by fitting stellar 
population synthesis models to the {Ha portion of the 
sources') optical spectra. Combining the stellar ages with 
2MASS K-band luminosities (where available) we com- 
puted the stellar masses of our sources following Drory et 
al. (2004). Drory et al. have parameterized the mass-to- 
light ratio in K-band as a function of stellar age (see their 
Fig. 1) using simple stellar population models (Maraston 
1998), and a Salpeter initial mass function. The total 
systematic uncertainty of such a derived mass-to-light 
ratio is estimated to be ^ 25 — 30%. Lastly, from the 
C0(1— )>0) luminosities inferred for our sources (see tabxo 
) we estimated the molecular {H2) mass using a conver- 
sion factor of a = M-h^/L'^^q — 1.5 (K km s^"^ pc^)"-'^ 
(|Evans et al.|[2005l) . 

We find systematic differences in the average black hole 
and host galaxy properties of the low- and high- exci- 
tation sources (i.e. LINERs and Seyferts, respectively) 
in our sample (tab:physprops ). This is illustrated in 
Fig. 131 where we also indicate the average properties 
of our low- and high-excitation radio AGN, computed 
using the ASURV statistical package and assuming log- 
normal distributions in luminosity and mass. The aver- 
age properties are specifically given in tab:averageprops 
□ Compared to LERAGN, HERAGN on average 

^ It should be kept in mind that in ASURV there is an imphcit 
assumption that the censored data follow a similar distribution to 
that of the measured population. If this is not the case, "average" 
values calculated by ASURV will be generally biased upwards (as 
our upper limits typically lie towards the bottom-end of the dis- 
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TABLE 1 

Physical properties of the z < 0.1 sample of radio AGN 



name 


rcdshift 


type 


il78-MHz 


-^'0.5-10kGv/-^^EDD 


stellar age 


M. 


Mbh 


Mh2 








[W/Hz/srad] 




[Gyr] 


[Mq] 


[Mq] 


[Mq] 


3C 31 


0.017 


Seyfert 


9.08x10^^ 


< 2.0 X 10^" 


3 


2.4 X 10** 


7.8 X 10^ 


(5.1 ±0.4) X 10** 


3C 33 


0.060 


Seyfert 


3.95 X 10^5 


1.6 X 10-^ 


5 


1.3 X 10** 


4.8 X 10** 


(3.75 ± 1.5) X 10** 


3C 98 


0.030 


Seyfert 


8.75x102-1 


3.7 X 10--* 


2 


7.9 X 10*0 


1.7 X 10** 


< 7.8 X lO'^ 


3C 321 


0.097 


Seyfert 


2.6x10^5 




13 


7.0 X 10** 




(3.3 ±0.6) X 10^ 


3C 403 


0.059 


Seyfert 


3.5x10^5 


3.3 X 10-3 


5 


2.4 X 10** 


2.6 X 10** 


(6.6 ± 1.6) X 10^ 


3C 449 


0.017 


Seyfert 


6.51x10^3 


< 7.0 X 10-3 


3 


2.4 X 10*0 


5.1 X 10^ 


(1.1 ±0.2) X 10* 


3C 452 


0.081 


Seyfert 


7.54x10^^ 


3.3 X 10-3 


13 


4.5 X 10** 


3.5 X 10** 


8.1 X lO**''^ 


Cen A 


0.0008 


Sey 2" 


5.4x10^3 


3.0 X 10-5 


- 


- 


2.0 X 10** 


1.4 X 10** 


3C 405 


0.0565 


Sey 2° 


4.90x10^7 


8.5 X 10-4 






2.5 X 10'-> 


< 3.3 X 10** 


3C 66B 


0.022 


LINER 


2.21x1024 


< 4.4 X 10-5 


3 


3.0 X 10*" 


6.9 X 10** 


< 7.8 X lO'^ 


3C 84 


0.018 


LINER 


3.74x102'' 


< 9.2 X 10-® 






1.9 X 10^ 


(2.14 ±0.02) X 10^ 


3C 264 


0.022 


LINER 


2.20x102'' 


< 1.8 X 10-5 


13 


4.4 X 10** 


7.1 X 10* 


(9.3 ± 1.8) X 10^ 


3C 272.1 


0.004 


LINER 


3.1xl022 


< 8.5 X 10-^ 


13 


2.9 X 10** 


1.5 X 10^ 


(9.3 ± 3.2) X 105-'= 


3C 274 


0.004 


LINER 


3.4xl024 


< 4.3 X 10-^ 


2 


1.5 X 10** 


2.4 X 10^ 


(1.65 ±0.15) X 10^ 


3C 296 


0.025 


LINER 


1.43x1024 


< 1.2 X 10-5 


13 


1.1 X 10*2 


1.3 X 10^ 


< 5.7 X lO'^ 


3C 338 


0.032 


LINER 


8.63x102-* 


< 2.0 X 10-5 


13 


1.3 X 10*2 


1.7 X 10'^ 


3 X 10^ 


3C 388 


0.091 


LINER 


4.29xl02-"5 
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< 1.2 X 10'-*''' 


3C 465 


0.030 


LINER 


6.41x102-* 


< 2.2 X 10--* 


13 


1.0 X 10*2 


2.1 X 10'^ 


< 1.95 X 10'' 


NGC 1265 


0.027 


LERG" 


3.39x102-* 


< 6.8 X 10-® 


13 


1.5 X 10*° 


1.0 X 10'^ 


< 5.7 X lO'^ 


NGC 6109 


0.0296 


LERG" 


1.86x102'* 










(1.3 ±0.3) X 10** 


NGC 6251 


0.0244 


LERG'' 


1.20x102-* 


< 2.0 X 10--* 






6.0 X 10** 


< 7.5 X lO'^ 



" Based on NED; LERG abbreviates "low-excitation radio galaxy" 

^ Adopted from Saripalli et al. (2007), and scaled to the cosmology used here. 

Not considered in our statistical analysis (see tabxo for details). 
The first, second and third columns denote the source, its redshift and AGN type. The last was inferred either via optical diagnostic diagrams 
(see Fig.[lJ, or adopted from N ED. The fourt h column shows the radio continuum luminosity, adopted from Evans et al. (2006). The fifth 
column, also adopted from Eva ns et al.l II2006I ). represents the accretion efficiency (in Eddington units) derived from X-ray observations of 
the cores of the AGN (the upper limits are obtained assuming Nh = 102-* atoms cm 2j see Evans et al. 2006 for details). The sixth column 
shows the stellar age of the source based on fitting stellar population synthesis models to the optical spectra of the sources (encompassing 
the Ha portion of the spectrum; see Tab. 5 in Buttiglione et al. 2009). The seventh column shows the stellar mass derived using the 
2MASS K-band luminosity and stellar age (where available) following Drory et al. (2004; see text for details). The second to last column 
shows the black hole mass, adopted from Evans et al. (2006), and the last column reports the molecular gas mass obtained from CO(l— >0) 
observations (see tabxo ) using a conversion factor of a = Mu^l ^'co ~ ^'** pc^)~^, and assuming Hq = 70, Om = 0.3, 

SIa = 0.7. The horizontal lines separate Seyferts (i.e. HERAGN; top) and LINERs (i.e. LERAGN; bottom). 



TABLE 2 

Summary of observations with the CARMA interferometer. 



Source 


RA (J2000) 


DEC (J2000) 


configuration 


obs. frcq. 


On-source 


beam 


rms / channel 










[GHz] 


time [hr] 


[arcsec] 


[mjy] 


3C 296 


14 16 52.94 


+10 48 26.50 


D & E 


112.460 


14.1 


6.4" X 5.5" 


2.1 


3C 321 


15 31 43.45 


+24 04 19.10 


E 


105.079 


4.0 


9.8" X 7.2" 


3.8 


3C 33 


01 08 52.86 


+13 20 13.80 


D & E 


108.746 


14.8 


8.7" X 6.3" 


1.8 


3C 403 


19 52 15.80 


+02 30 24.47 


E 


108.849 


5.7 


9.0" X 6.9" 


4.7 


3C 452 


22 45 48.77 


+39 41 15.70 


D & E 


106.634 


15.6 


6.8" X 5.1" 


1.8 


3C 465 


23 38 29.52 


+27 01 55.90 


E 


111.914 


8.9 


9.3" X 6.5" 


3.0 


3C 66B 


02 23 11.41 


+42 59 31.38 


E 


112.790 


4.4 


8.7" X 6.3" 


5.1 


3C 98 


03 58 54.43 


+10 26 03.00 


E 


111.914 


4.8 


9.3" X 6.7" 


3.7 


3C 83.1B (NGC 1265) 


03 18 15.86 


+41 51 27.80 


D & E 


112.241 


13.5 


4.5" X 3.7" 


2.1 


NGC 6109 


16 17 40.54 


+35 00 15.10 


D & E 


111.957 


10.6 


5.0" X 4.0" 


2.7 


NGC 6251 


16 32 31.97 


+82 32 16.40 


E 


112.526 


6.53 


9.0" X 8.2" 


3.5 
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Fig. 1. — Spectroscopic diagnostic diagrams for our 3CRR (filled dots) sources that separate AGN into LINERs and Scyfcrt s. The 
emission line fluxes have been taken from Buttiglione et al. (2009; see their Tab. 1). The regions (separated following Ka u ffmann et al.l 
[2003; Kcwlcv ct al. 20 0fl . [20061 1 that are occupied by star forming (SF), composite (Comp.) and AGN (Seyfert and LINER) galaxies are 
indicated in the panels. Each 3CRR galaxy is also labeled. 



TABLE 3 
Molecular gas properties 



Source 


z 


'S'cont 

[mjy] 


ZCO 


AtJpwHM 

[ km s-i] 


-fcO(i->0) 
[Jy km s"-'] 


[K km s"-"- pc^] 


3C 33 


0.060* 


31.8 ± 0.4 


0.060 


400 ± 200 


1.5 ±0.6 


(2.5 ±1.0) X 10* 


3C 66B 


0.022* 


113.6 ±0.8 






< 2.4 


< 5.2 X 10^ 


3C 83.1B (NGC 1265) 


0.027* 


40.5 ±0.4 






< 1.2 


< 3.8 X 10^ 


3C 98 


0.030* 


8.0 ±0.4 






< 1.3 


< 5.2 X 10^ 


3C 296 


0.025* 


144.4 ±0.4 






< 1.4 


< 3.8 X 10^ 


3C 321 


0.097* 


9.1 ± 0.6 


0.097 


320 ± 70 


5.0 ±0.9 


(2.2 ±0.4) X 10^ 


3C 403 


0.059* 


4.6 ±0.5 


0.059 


350 ± 100 


2.8 ±0.7 


(4.4 ±1.1) X 10* 


3C 452+ 


0.081* 


44.7 ±0.6 








< 5.4 X 10* 


3C 465 


0.030* 


20.2 ± 0.1 






< 0.3 


< 1.3 X 10^ 


NGC 6109 


0.0296** 


17.7 ± 0.3 


0.0301 


230 ± 50 


2.2 ±0.4 


(8.8 ± 1.7) X 10^ 


NGC 6251 


0.0244** 


624.9 ±0.6 






< 1.9 


< 5.0 X 10^ 


Cen A'^ 


0.0008** 










9.4 X lO'^ 


3G 272.1 (M84) 


0.004* 




0.0028 


200 


1.8 ±0.6 


(6.2 ±2.1) X 10^ 


3C 274 (MS?)* 


0.004* 




0.0035 


200 


20 ±2 


(1.1 ±0.1) X 10^ 


3C 31" 


0.017* 




0.0169 


450 


27 ±2 


(3.4 ±0.2) X 10* 


3C 449*" 


0.017* 




0.0169 


500 


6± 1 


(7.6 ± 1.3) X 10^ 


3C 84 (NGC 1275)" 


0.018* 




0.0176 


200 


104 ± 1 


(1.43 ±0.01) X Iff 


3C 264* 


0.022* 




0.02 


200 


3.5 ±0.7 


(6.2 ± 1.2) X 10^ 


3C 338'* 


0.032* 




0.030 




0.495 


2.0 X lO'^ 


3C 405" 


0.0565** 








< 1.5 


< 2.2 X 10* 


3C 388 


0.091 













The columns show the source, its redshift, the observed continuum flux density (S'cont), the redshift based on the CO(l— >-0) emission line 
(^co)i the line width at half-maximum (i^fwhm)! the CO line intensity (JcO{l->0))) ^^'^ luminosity {L'qq\ see eq. 4 in Evans et al. 2005). 
For sources in which the CO line was not detected we report 3(t upper limits, computed assuming AiipwHM = 300 km s~^. 
* adopted from Buttiglione et al. (2009) 
** adopted from Evans et al. (2006) 

+ Due to the strong contribution of complex, steeply sloped mm continuum emission from extended radio jets to the mm emission of this 
source, the continuum was fitted over only 33 channels, where the jet contribution is estimated to be small after deconvolution. Due to 
this uncertainty, however, we do not consider this source in our statistical analysis. 

" 2coi AtipvvHM and /cO(i— >0) adopted from Evans et al. 2005. L'qq was computed using the cosmology adopted here. 
^ zcOi Ai)pwHM and /co(i— >0) are taken from Flaquer et al. (2010). Given that their observations were conducted with the IRAM 30 m 
telescope, we take 1 K = 4.95 Jy, and compute L'qq using the cosmology adopted here, 
tentative detection 

7qo(i— >0) adopted from Leon et al. (2001). 
" adopted from Eckart et al. (1990), and scaled to the cosmology used here. 
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Fig. 2. — CO(l— >0) spectra (histograms) of the 11 radio galaxies in our sample observed with CARMA. The data are shown at a resolution 
of 31.25 MHz (83-89 km s"!). The black curves arc Gaussian fits to the line profiles and rest-frame 2.6 mm continuum emission. 
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TABLE 4 

Average properties of the z < 0.1 high- and low- excitation radio AGN 



AGN 


-f^l78-MHz 


stellar age 


M* 






type 


[W/Hz/srad] 


[Gyr] 


[Me] 


[Mq] 


[Mq] 


HEHAGN 


(7.2 i 1.9) X 10-4 


6. .3 ± 1.6 


(1.7 ±0.7) X 10** 


(2.-5 ± 1.0) X 10** 


(2.9 ± 1.2) X l(f 


LEliACiN 


(2.. J = 1.1) X 10--* 


10.2 - 1.1 


(2.1± l.i) X 10** 


(1..3 - 0.2) X 10'' 


(!.:-! - 1.9) X lO'-' 



* The given limit was computed excluding the tentative CO detection in 3G 272.1 (see tabxo ). Including the gas mass for this source 
yields an average of (1.8 ± 1.5) x lO'' Mq. 
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Fig. 3. — Distribution of physical properties of high- and low-excitation radio AGN (HEARGN and LERAGN, resp.), drawn from 
tab:physprops (excluding the tentative CO detection in 3C 272.1). Average values (given in tab:averageprops ) are indicated by vertical 
dot-dashed lines 
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have a factor of ~ 3 higher radio continuum luminosi- 
ties, significantly higher accretion efficiencies, but about 
an order of magnitude lower mass central black holes. 
Furthermore, their host galaxies have about a factor of 
1.5 younger stellar populations and stellar masses, but 
about a factor of ~ 7 higher molecular gas masses. As 
discussed in the next section, this is consistent with the 
idea that high- and low- excitation radio AGN form two 
physically distinct populations of galaxies that reflect dif- 
ferent phases of massive galaxy formation. 

4. DISCUSSION AND SUMMARY 

Our main result is that HERAGN have systemati- 
cally higher molecular gas masses (a factor of ~ 7; see 
tab:averageprops ), compared to LERAGN. Flaquer et 
al. (2010) have found a similar trend by dividing their 
sample (~ 50 radio AGN observed with the IRAM 30m 
telescope, partially overlapping with our sample) into FR 
class I and II objects. They find that the molecular gas 
mass in FR lis is a factor of ~ 4 higher than that in 
FR Is. The FR class can be taken to roughly correspond 
to the low- and high- excitation classification^ Flaquer 
et al. (2010) have, however, concluded that the system- 
atic differences they find are likely a result of a Malmquist 
bias, i.e. simply due to a systematically higher redshift 
of their FR-II sources. Although our HERAGN lie on 
average at a slightly higher redshift, compared to our 
LERAGN (0.046 vs. 0.030, resp.) in the following we ar- 
gue that the systematic differences we find in molecular 
gas mass are not due to a Malmquist bias. 

Moric et al. (2010) have shown that the redshift dis- 
tributions of carefully selected samples of radio-selected 
LINERs and Seyferts are approximately the same (see 
their Fig. 6). This eliminates Malmquist bias from their 
results. They find that the detection fraction in the FIR 
is significantly lower for LINERs than for Seyferts (6.5% 
vs. 22%, resp.) in their sample. Assuming that the 
star formation law parameterized by L'qq (as a proxy 
for total gas mass) and Lyyr. (as a proxy for star forma- 
tion rate; e.g., Kennicutt 1998; Solomon & Vanden Bout 
2005; Bigiel et al. 2008), on average, correctly represents 
the star formation properties of these samples (as con- 
firmed by the CO /FIR luminosities of the IRAS detected 
sources analyzed here; see Fig. [4|), the lower average FIR 
luminosity in low-excitation sources (i.e. LINER) implies 
lower gas masses than in high-excitation (i.e. Seyfert) 
types of galaxies. A similar result is obtained based on 
average (optically derived) star formation ratefl, sug- 
gesting that those in LINERs are by about a factor of 
3 lower than in Seyferts in a redshift-matched sample. 
These findings suggest that the systematic differences 
in molecular gas mass in high- and low-excitation radio 
AGN are physical, and not due to Malmquist bias. 

The systematically higher molecular gas masses that 
we find in HERAGN, relative to LERAGN in our z < 0.1 

* Almost all FR I - low power - radio galaxies are LERAGN, 
while optical hosts of FR lis, which are typically more powerful 
than FR Is (Fanaroff & Riley 1974; Owen 1993; Ledlow & Owen 
1996), usually have strong emission lines. Note however that the 
correspondence between the FR class and the presence of emission 
lines is not one-to-one. 

^ Moric et al. (2010) derived star formation rates for each galaxy 
in their sample via stellar population synthesis model fitting to the 
SDSS photometry of the host galaxy (see also Smolcic et al. 2008). 
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Fig. 4. — CO vs. FIR luminosity for our local AGN sources 
detected with IRAS. The lines represent the L'qq—Lyyr, correlation 
derived by IHiechers et al.l 120061) . 

radio AGN sample, are in excellent agreement with the 
systematic differences in various properties of high- and 
low- excitation radio AGN, both on pc- and kpc- galaxy 
scales (see Sec.[T]and tab:physprops ). 

We find that, on average, HERAGN have lower stel- 
lar masses and stellar ages compared to LERAGN (see 
tab:averageprops ; see also Smolcic 2009). This is con- 
sistent with HERAGN and LERAGN being green val- 
ley and red sequence sources, respectively. Furthermore, 
we show that HERAGN have on average higher radio 
luminosities than LERAGN, consistent with the results 
presented in Kauffmann et al. (2008). Kauffmann et al. 
have shown that the fraction of radio AGN with strong 
emission lines in their spectra significantly rises beyond 
10^5 W Hz^^. In general, the comparison of the 
black hole and host galaxy properties inferred for our 
21 z < 0.1 AGN with much larger samples of radio AGN 
(Kauffmann et al. 2008; Smolcic 2009) suggests that our 
AGN sample is representative of high- and low-excitation 
radio AGN in the nearby universe. 

From the average stellar masses that we infer for our 
high- and low excitation sources we extrapolate that 
they occupy - 3 x 10" Mf^ and ~ 5 x 10^^ Mf^ ha- 
los, r espectively (e.g. IBehroozi et al.|[20Tot IMoster et al.l 
|2010| ). Compared to the systematic molecular gas mass 
difference, this yields an even more dramatic discrep- 
ancy of more than 2 orders of magnitude in the average 
molecular gas fractions in HE- (^ 10^^) and LE-RAGN 
(~ 9 X 10~^). The discrepancy remains significant (about 
an order of magnitude) if the average gas-to-stellar mass 
fraction (which can be interpreted as star formation effi- 
ciency) is considered. 

On small scales, the average black hole accretion effi- 
ciencies in HE- and LE-RAGN suggest different super- 
massive black-hole accretion mechanisms (standard disk 
accretion of cold gas in HERAGN vs. Bondi accretion 
of hot gas in LERAGN; see Evans et al. 2006; Hard- 
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castle et al. 2006). Furthermore, the higher black hole 
masses in LERAGN suggest a later evolution stage of 
their host galaxies, compared to that of HERAGN. This 
is further strengthened by the higher stellar masses in 
LERAGN, as well as older stellar ages, and less mas- 
sive gas reservoirs. In the blue-to-red galaxy formation 
picture, blue gas rich galaxies are thought to transform 
into read-and-dead gas-poor galaxies, the stellar popula- 
tions in the host galaxies of HERAGN are expected to 
be younger and have lower masses, while their molecular 
gas reservoirs - fueling further stellar mass growth - are 
expected to be higher than those in LERAGN. This is 
in very good agreement with the results presented here. 
Thus, in summary, our results strengthen the idea that 
low- and high-excitation radio AGN form two physically 
distinct galaxy populations that reflect different stages 
of massive galaxy formation. 

The authors thank F. Bertoldi and K. Knudsen for in- 



sightful discussions. The research leading to these results 
has received funding from the European Union's Seventh 
Framework programme under grant agreement 229517. 
DR acknowledges support from NASA through an award 
issued by JPL/Caltech, and from NASA through Hub- 
ble Fellowship grant HST-HF-51235.01 awarded by the 
Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astron- 
omy, Inc., for NASA, under contract NAS 5-26555. Sup- 
port for CARMA construction was derived from the Gor- 
don and Betty Moore Foundation, the Kenneth T. and 
Eileen L. Norris Foundation, the James S. McDonnell 
Foundation, the Associates of the California Institute of 
Technology, the University of Chicago, the states of Cal- 
ifornia, Illinois, and Maryland, and the National Science 
Foundation. Ongoing CARMA development and oper- 
ations are supported by the National Science Founda- 
tion under a cooperative agreement, and by the CARMA 
partner universities. 



REFERENCES 



Allen, S. W., Dunn, R. J. H., Fabian, A. C, Taylor, G. B., & 

Reynolds, C. S. 2006, MNRAS, 372, 21 
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Behroozi, P. S., Conroy, C, & Wechsler, R. H. 2010, ApJ, 717, 

379 

Bell, E. F., et al. 2004a, ApJ, 608, 752 
Bell, E. F., et al. 2004b, ApJ, 600, Lll 
Borch, A., et al. 2006, A&A, 453, 869 

Bower, R. G., Benson, A. J., Malbon, R., Helly, J. C., Frenk, 
C. S., Baugh, C. M., Cole, S., Lacey C. G. 2006, MNRAS, 370, 
645 

Brown, M. J. I., Dey, A., Jannuzi, B. T., Brand, K., Benson, 
A. J., Brodwin, M., Croton, D. J., Eisenhardt, P. R. 2007, ApJ, 
654, 858 

Buttiglione, S., Capetti, A., Celotti, A., Axon, D. J., Chiaberge, 
M., Duccio Macchetto, F., & Sparks, W. B. 2009, 
arXiv:0901.1764 

Buttiglione, S., Capetti, A., Celotti, A., Axon, D. J., Chiaberge, 

M., Macchetto, F. D., & Sparks, W. B. 2010, A&A, 509, A6 
Buttiglione, S. et al., 2011, A&A, 525, A28 
Croton, D. J., et al. 2006, MNRAS, 365, 11 
Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615 
Drory, N., Bender, R., Feulner, G., Hopp, U., Maraston, C, 

Snigula, J., & Hill, G. J. 2004, ApJ, 608, 742 
Eckart, A., Cameron, M., Rothermel, H., Wild, W., Zinnecker, H., 

Rydbeck, G., Olberg, M., & Wiklind, T. 1990, ApJ, 363, 451 
Evans, A. S., Mazzarella, J. M., Surace, J. A., Frayer, D. T., 

Iwasawa, K., & Sanders, D. B. 2005, ApJS, 159, 197 
Evans, D. A., Worrall, D. M., Hardcastle, M. J., Kraft, R. P., & 

Birkinshaw, M. 2006, ApJ, 642, 96 
Faber, S. M., et al. 2007, ApJ, 665, 265 



Fanidakis, N., Baugh, C. M., Benson, A. J., Bower, R. G., Cole, 
S., Done, C, & Frenk, C. S. 2010, MNRAS, 1547 

Ocana Flaquer, B., Leon, S., Combes, F., & Lim, J. 2010, A&A, 
518, A9 

Hardcastle, M. J., Evans, D. A., & Croston, J. H. 2006, MNRAS, 
370, 1893 

Jackson, N., & Rawlings, S. 1997, MNRAS, 286, 241 
Kauffmann, G. et al., 2003, MNRAS, 341, 33 

Kauffmann, G., Heckman, T. M., Best, P. N. 2008, MNRAS, 384, 
953 

Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., 

Trevena, J. 2001, ApJ, 556, 121 
Kewley, L. J., Groves, B., Kauffmann, G., Heckman, T., 2006, 

MNRAS, 372, 961 
Leon, S., Lim, J., Combes, F., & van-Trung, D. 2001, QSO Hosts 

and Their Environments, 185 
Moric, I., Smolcic, V., Kimball, A., Riechers, D. A., Ivezic, Z., & 

Scoville, N. 2010, ApJ, 724, 779 
Moster, B. P., Somerville, R. S., Maulbctsch, C, van den Bosch, 

F. C, Maccio, A. V., Naab, T., & Oser, L. 2010, ApJ, 710, 903 
Peacock, J. A. 1985, MNRAS, 217, 601 
Riechers, D. A., et al. 2006, ApJ, 650, 604 
Saripalli, L., & Mack, K.-H. 2007, MNRAS, 376, 1385 
Sijacki, D., & Springel, V. 2006, MNRAS, 366, 397 
Sijacki, D., Springel, V., di Matteo, T., & Hernquist, L. 2007, 

MNRAS, 380, 877 
Smolcic, v., et al. 2008, ApJS, 177, 14 
Smolcic, V. 2009, ApJ, 699, L43 
Smolcic, v., et al. 2009, ApJ, 696, 24 

Solomon P. M., Vanden Bout P. A., 2005, ARA&A, 43, 677 
Willott, C. J., Rawlings, S., Blundell, K. M., Lacy M., & Eales, 
S. A. 2001, MNRAS, 322, 536 



